Irrotational dust with div H = 

W. M. Lesame 1 ' 2 , G. F. R. Ellis 1 , and P. K. S. Dunsby 

1 Department of Applied Mathematics, University of Cape Town, 

Rondebosch 7700, Cape Town, South Africa. 

2 Department of Applied Mathematics, University of Fort Hare, 

Private Bag XI 314, Alice 5700, South Africa. 

February 7, 2008 
Abstract 

For irrotational dust the shear tensor is consistently diagonalizable 
with its covariant time derivative: a a b = = & a b, a ^ b, if and only 
if the divergence of the magnetic part of the Weyl tensor vanishes: 
div H = 0. We show here that in that case, the consistency of the 
Ricci constraints requires that the magnetic part of the Weyl tensor 
itself vanishes: H a b = 0. 

Subject headings: 
cosmology - galaxies: clustering, 
formation - hydrodynamics - relativity - exact solutions 

1 Introduction 

The local non- linear dynamics of irrotational dust with a purely "electric" 
Weyl tensor has been investigated by Matarresse et al. [7], |[ ||. These 
authors characterized such solutions as 'silent universes', because the time- 
evolution equations form a set of ordinary differential equations. In it 
was shown that all the integrability conditions for this case are satisfied 
and are consistent, provided the initial constraint equations are satisfied. 
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Then the development along each world line is independent of those along 
neighbouring world lines; once the evolution is underway, no information is 
exchanged between neighbouring particles. 

The dynamics of such universes have been studied mostly by employing 
an orthonormal tetrad, which is a simultaneous eigentetrad of the shear a a b 
and the "electric" part of the Weyl tensor E^. Variables defined in terms of 
these tetrads allow development of useful phase -plane portraits @, which 
show inter alia that the generic collapse configuration is a cigar or spindle 
rather than a pancake. 

In the linear theory, the magnetic part of the Weyl tensor H a b only con- 
tains vector and tensor modes |]]; also if vorticity is zero, no vector modes are 
present, and H a b at most contains tensor modes with vanishing divergence, 
which describe gravitational waves. In the non- linear theory the physical 
content of H a b is however less clear. At the core of the dynamics of silent 
universe is the assumption that in the absence of H a b, no gravitational waves 
occur. 

For irrotational dust the shear tensor is consistently diagonalizable with 
its covariant time derivative: o~ a b = = cr ab , a ^ b, if and only if the 
divergence of the magnetic Weyl tensor vanishes. Furthermore if the electric 
Weyl tensor vanishes then the magnetic Weyl tensor vanishes ||. So for 
irrotational dust there are no solutions with a purely magnetic Weyl tensor. 

The aim of this paper is to prove that for irrotational dust there are no 
consistent solutions with a non-vanishing magnetic part of the Weyl tensor 
that has a vanishing divergence. Thus the only consistent solutions with 
vanishing divergence of H ab are those described by the silent universes [2,7- 
9]. Mathematically, we start from the requirement that the shear tensor o~ a b 
is consistently diagonalizable with & a b. We then prove that this is possible 
iff divH = (0). This is a dynamical restriction, for it is a constraint 
on the gravitational tidal field. The implication of this is seen in the time 
derivative of the Ricci constraints (which prescribes spatial restrictions on 
both variables a a b, H ab ), and this yields H ab = 0. 

Notation: Latin indices run from to 3, and Greek indices from 1 to 3; 
semicolons denote covariant derivatives. Covariant differentiation along the 
velocity vector u a is denoted by ( ) e.g., the acceleration vector is a a = ii a := 
u a -bU b . Kinematic and dynamic quantities are the same as in the other papers 
of this series (they are comprehensively defined in Q). 
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2 Equations for Kinematic and Dynamic 
Quantities 



2.1 Conservation Equations 

For a perfect fluid, the fluid acceleration is only determined by pressure 
gradients so the restriction of vanishing pressure implies that a a = 0. This 
means that each fluid element moves along a geodesic. The conservation of 
energy and momentum T ab -b = leads to only one equation, the continuity 
equation: 

p = - P e. (f ) 

2.2 The Ricci identities 

With the second restriction of vanishing vorticity u a b = 0, the equations 
for the kinematic quantities follow from the Ricci identity: u a . d . c — u a]C] d = 
R a bcdU b - 

2.2.1 Propagation equations 

The expansion scalar G obeys the Raychaudhuri equation: 

+ |0 2 + 2a 2 + \k P = , (2) 

where a 2 = \(? ab <3 a b is the shear scalar and k = 8nG is the gravitational 
constant. The remaining kinematic evolution equation is for the shear and 
is given by: 

°ab + v ac v c b - l^ 2 h ab + |6a afe + E ah = , (3) 

where E ac = Ei ac \ = C a i> c dU b u d is the "electric" part of the Weyl tensor 
Cabcd (satisfying E ac u c = 0, E a a = 0). E ab is that part of the gravitational 
field which describes tidal interactions. The Weyl tensor can be decomposed 
into E a b and another tensor called the "magnetic" part: H ac = if( ac ) = 
\Vab 9h C ghcd u b u d (satisfying H ac u c = 0, H a a = 0). This is the part of the 
gravitational field that describes gravo- magnetic effects, and gravitational 
waves. 
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2.2.2 Constraint equations 



Besides the evolution equations for the kinematical equations, there are sev- 
eral constraints that our variables must satisfy. On setting p = uo ab = we 
obtain as non - trivial constraints, 

h\ (f Q' b - h d c a bc , d ) = , (4) 
the '(0, z/)' field equations, and 

H ad = -rfah'daff'^fbcu* ■ (5) 

which we refer to as the l H ab constraint, also (^) and (|5|) are both referred 
to as Ricci constraints. As we are considering the zero - vorticity case, we can 
also write down the Gauss- Codacci equations for the 3 -spaces orthogonal 
to u a (see (U); however we do not need them for what follows. 

2.3 Bianchi identities 

Additionally, the Bianchi identities must be satisfied, as they are the inte- 
grability conditions for the other equations. 

With the restrictions we have put on so far, they take the form: 

h\h d s E as . d - v tbp i Ub <T d q H dq = \h\p< b , (6) 

h m a h t c E ac + h mt a ab E ab + QE mt + J mt - 3E s {rn a t)s = -\pa tm , (7) 
which we refer to as the "div E" and equations respectively, and 

h a H as . d h d s +V tbpq u b (T d p E qd = 0, (8) 

/i m a /i* c ir c + h mt a ab H ab + QH mt - I mt - 3H s im (T t)s = , (9) 

which are the "div H" and "if" equations respectively, where we have defined 
the curls of E and H respectively as 

'Curl E' : I mt = -h a {m ri t)rsd u r E a S]d (10) 

'Curl H' : J mt = h a {m rf )rsd u r H a s . d . (11) 

The Bianchi identities are analogous to Maxwell's electromagnetic equation 
[f|]. The gradient \h b a p^ acts as a source of the divergence of the E ab field 
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in the "div E" constraint. For zero vorticity the H ab field is source -free in 
the "div H" constraint. In the case of the time - development equations the 
"curl E" term I mt acts as a source of H. On the other hand \pa tm acts as a 
source of E, as well as the "curl H" source term J mt . 



3 Tetrad Description 

We now show that the shear tensor a a b is diagonalizable in the same principal 
frame as a a b iff div H a t, = 0. 

We first introduce an orthonormal tetrad that diagonalizes the shear ten- 
sor i.e, 

a ab = (a^b). (12) 

This immediately implies da a b/dr = for (a ^ b), where u = d/dr, but this 
gives no direct restrictions on & ab because of the rotation coefficient terms 
in those quantities. However from (|3|) the tetrad satisfying ( [P^) also satisfies 
&ab = for a 7^ b provided 

E ab = (a^b) (13) 



in this tetrad. From the "div H" equation (§), equations (|i~2"D and fll3|) imply 
that 

dtvH = h t a H as . d h d s = 0. (14) 

Conversely if div H = then a a b and E a b are simultaneously diagonalizable 
from (pp, as shown in 0; furthermore then a a b = for a ^ b from (|3|). 
Hence, For irrotational dust the shear tensor is consistently diagonalizable 
with its covariant time derivative: a a b = = <j a & for a ^ b, if and only if the 
divergence of the magnetic Weyl tensor vanishes. 

We now adopt the assumption that div H = 0, as defined by equation 
(|I4"D; and choose the preferred frame that diagonalizes a a b and a a b- For 



simplicity we denote a aa (no sum) by a a , and E aa by E a . Both the shear 
tensor a ab and the 'electric' part of the Weyl tensor E ab satisfy the trace-free 
property, i.e., 

Ei + E2 + -£"3 = , 
ci + 02 + 0-3 = . (15) 
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If we now write out = & ab = cr ab - c u c for (a ^ b), we get the following tetrad 
relation: 

(<r - <7 6 )r°o6 = 0; a^b (16) 

where u c = 5 c q and u a u a — — 1. For an arbitrary shear tensor one may deduce 
from (|16|) that 

r a 06 = (a^b), (17) 

and hence the tetrad is Fermi-propagated along u a . Condition ([17]) can also 
be shown to be valid for the case of degenerate shear. This is achieved by 
performing a rotation in the degenerate plane and using this tetrad freedom 
to obtain ([[7]), see [6]. Furthermore for dust the vanishing of vorticity in a 
shear eigenframe is equivalent to the conditions 

T° ab = . (18) 

Now by (|17D and the diagonality of E ab , it follows that also E ab is diagonal: 

E ab = {E a - E b )T a ob = 0, a + b. 

Then from the time development equation (|7|) it follows that "curl H" is also 
diagonal i.e., 

jmt = hJ m v t)rsd u r H a s . )d = , m^t. (19) 

Consistent with the rest of the notation, we denote J aa by J a . Then the fact 
it is trace free becomes Ji + J 2 + J 3 = 0. 



3.1 Propagation equations 

Equations (|l])-(0) governing the evolution of irrotational dust may now be 
written as a set of propagation equations tied to a set of constraint equations. 

A direct conversion of equations ([!]), (Q), (|3|) and (|7|) into a tetrad system 
which is an eigenframe for both the shear tensor and the electric Weyl tensor 
yields the following time-evolution equations: 

P = -p& , 

= _l e 2_( (Tl 2 + (T22 + (73 2 ) _l p) 
<V = -K) 2 -|©^ + |((7i 2 + ( T 2 2 + ( T3 2 )- J E M , 

E^ = -6£ M - \po^ - J M + 3a M £ M - {o x E x + a 2 E 2 + cx 3 E 3 ) , (20) 
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where \i = 1,2,3 and no sum is carried out. The propagation equation for 
the magnetic Weyl tensor H a b follows from the Bianchi identity @ and has 
the form 

H n = -OH u + 3aii7 n + J u - ((Tii^n + a 2 H 22 + 0-3^33) , 

H22 = —QH22 + 3(T2-f^22 + ^22 — (Cl-f^ll + 0"2-f^22 + 0~ 3 H 33 ) , 

H 33 = -QH 33 + 3a 3 H 33 + I 33 - {atH n + a 2 H 22 + a 3 H 33 ) , (21) 

-^23 = —6-^23 — §0"l-£f23 + -?23 , 
H12 = —QH12 — ^<J 3 Hi2 + I12 , 

H 31 = -QH 31 -la 2 H 31 + I 31 . (22) 

It is important to note that in each of these equations, the left hand side 
is in terms of the covariant derivative T = T. a u a , and so in general involves 
rotation coefficients; for example, in the obvious notation, the left-hand side 
of the 3rd equation fl2"0|) is 



o~ v = °v,o — 2cr I/ r I/ o^ {no sum) (23) 
and that of the 4th equation is 

E v = E vfl - 2E V T\ U {no sum). (24) 

3.2 Ricci Constraints 

The tetrad form of the Ricci constraints (El), fl5T) follows: the "(0,i/)" 

field equations (HI) take the form 



1^6 = <9i<7i + (<n - o- 2 )T 2 21 + {ai - a 3 )T 3 31 
|a 2 9 = d 2 a 2 + {cr 2 - o-^T 1 ^ + {0-2 ~ or 3 )T 3 32 



§<9 3 6 = d 3 a 3 + {a 3 -a 1 )T\ 3 + {a 3 -o-2)T 2 23. (25) 
The ll H a b" equations (|5|) take the form 



(26) 





= r\ 3 i 


>3 




- r 1 32 ((T2 




H22 


= r 2 3ll 


>1 


- ^2) 


- T 2 13(0"3 


- ^2) 


H 33 


= r 3 12l 


v°"2 


- 0"3) 


-r 3 21 (^ 


- 0-3) 
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#23 = \ [#1 (CT2 - O-3) + r 3 3 i((Ti - O3) - r 2 2l((Tl - <T 2 )] , 

#31 = I [<9 2 (<x 3 - oi ) + r 3 32 ((T3 - «r 2 ) - r 1 i 2 (^i - <T 2 )] , 

#12 = H 5 3K-02) + r 1 1 3(^i-^)-r 2 23 (a 2 - ( 73)] . (27) 



3.3 Form of H ab 

Given our choice of a tetrad frame which diagonalizes a ab (and hence accord- 
ing to theorem 1, divH = 0) we show below that consistency of the Ricci 
constraint requires that H ab = 0. 

To calculate time derivatives of the constraints we first start from the 
covariant form of the constraint equation, say (G) and then take the time 
derivative as defined by (G) = (G). a u a . We commute derivatives and use the 
constraint equations; the result is then converted to tetrad form. We give 
an example of such a calculation in Appendix A (using the ^H a b constraint). 
This approach avoids any direct calculations of the time derivatives of Ricci 
coefficients r a bc . 

The time derivative of the "(0, z/)" constraint (f|) in covariant form is 



= h 



eb 



(eg- 



h d c (a b c . 



If we use 

(eg- = 

and substituting into 



e 



Vb C ]P u P ;d + R c qp d (rb q u p - R q bpd a g c u p 
we get 



(28) 

(29) 
(30) 







h 



eb 



e) . -9/ ;l 



h d c {(a b c 



'■d 



h eb h d c {R c qpd a b q u p - R q bpd a q c u p } . 



(31) 



We may now convert (|3lD into tetrad form and use (p0|) and (p5|). The 
Riemann tensor is given in terms of the Weyl tensor in Appendix A. On 
further simplification (see that Appendix for details) we obtain 







{a 2 -<r 3 )H- 



23 







(0-3 -ai)H : 



31 







(01 - cr 2 ) H 



12 • 



(32) 
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The time derivative of the u H a b" constraint ([3j) in covariant form may be 
written as 



(33) 



H ad = h {t a H s \ v c fbs u f [(v b t). c - A-, P u p - c + R\ pc a\vP - R\ vc a l 

The tetrad form of (|33"D for a ^ b yields 

= aiH 23 , = a 2 H 31 , = a 3 H 12 (34) 
(see the Appendix B for details). 

For nonzero shear conditions (^) and fl3~4] ) yields 

H 12 = H 23 = H 31 = 0, (35) 

and thus the "magnetic" part of the Weyl tensor H a b is also diagonal. 
Through the H equations, this implies curl E too is diagonal. We henceforth 
write H aa as H a ; then the trace free property of H ab is 

H l + H 2 + H 3 = Q . (36) 

The tetrad form of (|33D for a = b yields (see Appendix B) 

0l H x = a 2 H 2 = a 3 H 3 , (37) 

We point out that diagonal property (^) of H a b was not used in obtaining 
(|37|). We now introduce a constant A that relates the shear eigenvalue a\ = o 
to <y 2 as follows 

o\ = a, o 2 = Act . (38) 

If the shear tensor is degenerate in the ei, e 2 plane then A = 1. Now equation 
( |37D o\H\ = a 2 H 2 prescribes the following relation on the eigenvalues Hi, H 2 
of the magnetic Weyl tensor: 

Hi = XH, H 2 = H . (39) 

The trace-free property yields 

<t 3 = -(1 + A)<7, H 3 = -(1 + X)H, (40) 
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so now from (|37| ) if we use <J\H\ = we obtain 

= (1 + A + \ 2 )<jH 



(41) 



The following cases satisfy fl4"T| ) 

Case A: H = with a 7^ and A remaining arbitrary. This case include 
degenerate shear, A = 1, and is studied in |7], |8], [| and |J. Variables defined 
in terms of the above tetrad frame allow development of useful phase - plane 
portraits, which show inter alia that the generic collapse configuration is a 
cigar or spindle rather than a pancake. 

Case B: o = 0. For this case all the shear eigenvalues vanish and this 
implies H = and E = and the model is FRW. 

Case C: (1 + A + A 2 ) = 0. The values of A are complex. The magnitudes 
of the shear tensor a mag and the magnetic Weyl tensor H mag have the tetrad 
form 

O-mag = \{o 2 + <T 2 2 + 0" 3 2 ). 

H mag = #! 2 + H 2 2 + H 3 2 . (42) 

On using (|38D and (|40|) for this case we get 

Vmag = (1 + A + A> 2 = 0, H mag = (1 + A + \ 2 )H 2 = (43) 

and hence both shear and the magnetic Weyl tensor are zero, leading to the 
same results B above. 

Thus we can formulate the following 

Theorem: 

For irrotational dust the divergence of the magnetic Weyl tensor 
vanishes, divH = (or equivalently the shear tensor is consis- 
tently diagonalizable) if and only if the magnetic Weyl tensor 
vanishes, H a b = 0. 

4 Conclusion 

For irrotational dust p ^ the existence of a tetrad frame which is a principal 
frame for the shear tensor a a b and its covariant time derivative a a b requires 
that the divergence of the magnetic Weyl tensor vanish. We have shown 
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here that if we employ this frame then the magnetic Weyl tensor itself van- 
ishes; H a b = 0, as a consistency requirement of the Ricci constraints. This 
establishes a new property, that the only consistent solutions for irrotational 
dust with vanishing div H are those described as the 'silent' universe @,[[|]. 
Hence gravitational waves interacting with irrotational dust will have to have 
div H ^ 0, contrary to the usual result of linearised theory. 

The key point here is that (assuming the fluid is irrotational), our result 
comes from the second term in the 'div H' equation (8); but when we linearise, 
that term is necessarily second order. Thus always DivH = (to first order) 
in the linear case. Hence our exact result comes from a term which plays 
no role in the linearised theory, if we discard all second order terms in all 
equations; in this case the surviving Riemann terms would be disregarded 
(as they always consist of a Weyl tensor component, which is first order, 
multiplying a first order quantity: see (32), (34), (37)). 

However this argument needs to be treated with care. One needs to recall 
that second order terms can only be dropped from an equation if there is a 
non-zero first order term present, so that the second-order term is negligible 
relative to the first-order one. However this argument cannot be applied 
to the key equation (8): both terms are second order, and as there is no 
first order term, we have to ensure that the equation is true to second order 
accuracy. Thus even in linearised theory, we cannot ignore this second-order 
equation. The following intriguing situation results: if we have a linearised 
solution where H a b ^ and div H is non-zero but second order, we can 
presumably get a consistent solution. If however H ab ^ with div H exactly 
zero, the solution will not be consistent - because of the above proof. 

Thus in the linear theory, it is possible to have models where Div H = 
to second order but H a b ^ 0, and indeed that is usually assumed for 
gravitational waves. Our result (for irrotational dust only) shows that in this 
case in fact we must have div H ^ 0, although it is second order. We have 
therefore an example of linearization instability in that the usual process of 
linearization leads to a different answer than the exact result - which in fact 
constrains the linearised solution, even though this is usually not commented 
on. 
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A Time propagation of the "(0, //)" constraint 

The covariant form of the "(0,/i)" constraint (4) is 




(44) 
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with tetrad form 



foe 



fa! + (ai - a 2 )r 2 21 + (o-! - (T 3 )r 3 



§d 2 e = d 2 a 2 + (o- 2 - (7 1 )r 1 12 ( ( j 2 - (7 3 )r 3 32 
§a 3 e = <v 3 + ((t 3 - ^or^a + - <r 2 )r : 



31 



23 ? 



(45) 



The covariant time propagation of ([44] ) is (31), 

= h eb [i(Q, b y -h d c (a b c . d y 

= h eb [(6) b - Q ;p u p . b - h d c {(a b c ), d - a b C ; P u p ., d } 
-h eb h d c {R c qpd a b q vP - R% d (T q c v?} . 



(46) 



We now write the Riemann tensor R s mpc in the terms of the Weyl tensor as 



R rape C m pc 2 (j? pRcm Q cRpm 9rapR c 9mcR p) 

' 9 c9pm) 



~jr{9 p 9c 



(47) 



where 



C 



mpc — (j] mij^lpckl ~\~ 9 mij9pckl)WU E 

+ (v S mij9pcki + rjpcki9 S mij)u l u k H hd 



(48) 



and 



9 mij — 9 i9mj 9 jQr. 



The two terms involving the Riemann tensor R a bcd in 
covariant form as 



(49) 

may be written in 



(rty 



h eb h d c [R c qpd a b q u p - R q bpd a q c u p 



-h eb h d c u p [\a\ (g c p R qd + g c d R qp - g qp R c d + g qd R c 

—\ aC q (g 9 P Rbd + g q dRbp — gb P R q d + 9bdR q p) 



R 



{g C p9qd - 9 C d9qp) + f ^ p9bd ~ ' d9bp) 



- h eb h d c u p u l u k E jl [a q b r] c q ijVpdki + & 9 b (g C i9qj ~ 9 C jg q i) (g P kgdi ~ 9pi9dk) 
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-v c q r) q bijV P dki - o c q {g q i9bj - g q jgbi) (g P kgdi - g P Wdk)} 



- h eh h d c u v u l u k H 31 [a\r] c qij (g pk g dl - g p ig dk ) ~ v c q V Q bij (g P kgdi - g P igdk) 
+(j q bV P dki (g c ig Q j - g c jg q i) - v c q v P dki (g q ig C j - g q jga)] ■ (50) 

The only non-vanishing terms in (|50f) are 

(RTf = -h eb h d c uWu k W l [a\ V c qij g pk g dl - a c qV \ i]9pk g dl } 

= rf qij a e mKv? - h eh rf hl3 a ql W l u l . (51) 

The tetrad form of fl5"TD becomes 

(RTf = rf m a e m\ - r]\ 0j a ql W l . (52) 

and from (|52|) for e = 1 we obtain 

(RT) 1 = rf q ^a x mK-rf m cy ql E^ 

= r 1 3 ia2 (J 1 H 2 3 + r ] 2 103 a 1 H 3 2 - v 3 w2^H 23 - V 2 W3 a 2 H 32 

= #23(03 - o- 2 ) • (53) 



Similarly from (|52|) for e = 2, 3 we get 



(RT) 2 
(RT) 3 



H 31 (a 1 - cr 3 ) , 
H 12 (cr 2 - cri) . 



(54) 



We now focus on the tetrad form of the time propagation equation 
for e = 1 ; in expanded form this becomes 







d x d - OidxO - [d x dx - 0i<9i<7i] 



~r 2 21 [(&X - CT 2 ) - 6 2 (<7l - Cr 2 )] - T 3 3i [(&x - & 3 ) - 3 (c7i - (J 3 )\ 

+H 23 (a 3 - a 2 ) , 



(55) 



where the last term on the right hand side is the contribution of the Riemann 
term (RT) 1 calculated in (|5^). On using and from equation (20) we get 







-§00i0 - dx (ax 2 + a 2 2 + a 3 2 ) - \d lP - (ax + \0)d x t 
-2axdxax - ~Qdx(Jx ~ \^id\0 + \dx(ox 2 + o 2 2 + a 3 2 ) - oxdxOx 
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r 2 

"J- 21 



1*1 



2 



-r 3 



31 



lg( ai - a2 ) - ( El - E 2 ) - (cx 2 + \e){a x - a 2 ) 
W - a 3 2 ) - U{a x - o 3 ) - {E x - E 3 ) - {a 3 + \9){jj x - a 3 ) 



+H 23 (a 3 - o 2 ) 



(56) 



where we also used u a . b = a ab + \6h a b and # M = o M + ^9. We now apply to 
(56) the "(0,/i)" constraint fl45|) and the Ll div E n constraint written here in 
tetrad form as 

\d lP = dtEi + (E 1 - E 2 )T 2 2l + (E y - E 3 )T 3 31 + H 23 {a 3 - a 2 ) , 
\3 2 p = d 2 E 2 + (E 2 ~E 1 )T 1 12 + (E 2 -E 3 )T 3 32 + H 31 (a 1 -a 3 ) , 
\d 3 p = 8 3 E 3 + (E 3 - E^T 1 ^ + (E 3 - E 2 )T 2 23 + H l2 {a 2 - a x ) 

to obtain the following form of 



(57) 







3<Ji9io-i — <9i((Ji 2 + o 2 + a 3 



(58) 



-r 2 2i(o"i - o 2 )(o 3 - a 2 ) - r 3 3 i(cr! - a 3 )(a 2 - a 3 ) . 
The first term on the right of (|58D simplifies as follows 
3cri<9icri - di(ai 2 + o 2 + a 2 ) 

= 3cri9i<Ti - 2a 1 dia 1 - <9i(o 2 2 + cr 3 2 ) 

= Oi^o-i - 9i(a 2 2 + o 3 2 ) 

= \d l {-a l 2 -2a 2 2 -2a 3 2 ) 

= \d 1 [{a 2 -a 3 ) 2 -2a 2 2 -2a 3 2 } 

= -\di(a 2 2 - 2o 2 a 3 + a 3 2 ) 

= -\di{a 2 - cr 3 ) 2 

= -(o 2 - o 3 )<9i(o 2 - o 3 ) , (59) 

so now (^J) in (|58"D gives 

= (02-03) 9i(cr 3 - o 2 ) + r 2 2 i(o-i - cr 2 ) - r 3 3 i(o-i - a 3 ) 

= (a 2 -a 3 )H 23 , (60) 

where in the last step we used the U H 23 " constraint (27) 

#23 = |[9i(o2-o 3 )-r 2 2i(o 1 -a 2 ) + r 3 3i(o 1 -a-3)] . (61) 
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Similar relations to fl6"0|) follows for e = 2, 3 in ( [46] ) these are 








(0-3 - CTi)#3i 

(0-1 - a 2 )H 12 



(62) 



Thus we have derived the required equations (32). 



B The time derivative of the constraint 

We first write the covariant form (5) of the u H ao n constraint as 

-h t a h s d cr b ( t ;c r] s)fbc u f . 



H, 



ad 



W a h S d U f 



CT b t;cVsfb C + CT b s - c 7] t fb 



If we use 

s;c) (^" s]\c ® s;p^ ;c ~t~ R avc@ R q svc& 

the time propagation of (|63"D becomes 



(63) 



(64) 



Had = -^ah^U 1 



//.,/,/ <j \a b t ) ]C - a h t , p u p ,c} + Vtfb c {{cr b s )-,c ~ v b s- )P u p - c } ] 



Vsfb 



{R h qpc a\vP - R% c a h q vP) + vt f b c {R h qp c° q s uv - R q spc a b q u p } } .(65) 



Using the Riemann tensor R s m pc expression as given in the previous Ap- 
pendix, the two terms involving the Riemann tensor R a bcd in (0) m ay be 
written in covariant form as 



(RT) 



ad 



W a h s d u f 



Vsfb 



{R b qpc a\u p - R\ pc a b qU v} 



+Vtfb c {R b qP cO q s u p - R q spc a b q u p } 



—\h t a h s 'du* f r] S f b c \o q t [g b p R qc + g b P R qp — g qp R b c + fiv-R fe P ) 



x -a b q {g\R tc + g\R tv - g tp R\ + g tc R\ 



o q t [g b P g q c - g b c g qP ) + f <A {g q P gtc - g q c gt P ) 



+ u l u k E jl {a q t r] b qij 7] pck i + o q t (g b ig qj - g b jg qi ) (g P kg c i - g P ig C k) 
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-V b qTi q tijTlpckl ~ (fl^tf ~ ^ j9ti) (9pk9cl ~ 9pl9ck)} 



+ u l u k E 31 {(r q t T] b qij (g P k9ci ~ 9 P i9ck) - (y b q if tij (g pk g c i - g p ig c k) 
+(^ q tv P cki (g b ig q j - g b jg q i) - (? b q ii P cki (g 9 igtj - g q jgu 



\h l a h s d u j 'r] t f b c \o q s (g b pR qc + g b p R qp — g qp R b c + g qc R b p 

2^ i (g q p-^sc g q cRsp gspR^c gscR q p) 
-%° q s (g b P g qc - g b c g q p) + f <A (g q p g sc - g q 



+ u l u k E jl {<j q sV b q ijV P cki + o q s (g b ig q j - g b jg q i) {g P kg c i - g P ig C k) 
' b q V 9 sijV P cki - cr b q (g q ig S j - g q jg S i) (g P kg c i - g P ig C k)} 



-a 



+ u l u k H 31 { a q s r] b qij (g pk g c i - g p ig c k) - cr b q 7] q sij (g pk g d - g P Wck) 

+v q sV P cki (g b tg q j - g b jg q i) - v b q v P cki (g q ig S j - g q jg sl )}] • (66) 

The only non-vanishing terms in (|66|) are those containing H-* 1 i.e, 

(RT)ad = -\h t a h s d ufr] sfb c u i u k H^ (o q t if qij g pk g d - a b q r] q tij g pk g cl ) 
-\h l a h s d u s ri tfh c u l u k H ]l [a q s r] b qij g pk g d - a b q r] q sij g pk g c ^j 
= -\h t a h s d u s u v u % u k H : > c r] sfb c (g pk r] b qij cr q t - g pk ri q ujO- b q ) 

+Vtfb c (g P kV l \ij a<l ' s - gpkV q sijO- b q ) ] 

= Itfah'jufuWc [ Vsfb c (rf qi] a q t - v q tijcr b q ) 

+r }tfb c u f {r ] b qij a\-r ] q sij a b q ) ] . (67) 



Equation (|67"D represent the covariant form of the Riemann term. The tetrad 
form of (p7[) for diagonal elements a = d becomes 



{RT)c 



VaOb C V b aOjH J c (a a - a b ) . 



and hence 



(RT) 



11 



Hsfa - cr 2 ) + H 2 (<Ti - 0-3) 



(68) 
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{RT) 2 2 = H x (a 2 - cr 3 ) + H 3 (a 2 - cr x ) , 
(RT) 33 = H 2 {a 3 - <7i) +#1(0-3 - (7 2 ) . 

The tetrad form of the time propagation equation (|65|) for a 



(69) 

d now becomes 



-VaOb 
+ (RT) aa 

becomes 



d c a\ - 9 c da b a + T b ca {(& a - & b )} - 9 c (a a - a h ) 



(no sum over) a . 



(70) 



For a = 1 fl7 
#ii = r 3 21 [((7i - & 3 



2 (oi " 0- 3 )] - T 2 3i [((7! - 02) - #3(01 - 02)] + (i£T) 



11 • 
(71) 



On using in (20) and constraint (26) we get 

H n = -eH n + I n + {RT) n . (72) 
If we compare with the propagation equation for H xx in (21) i.e., 

H u = -6H 1X + I n + 3a x H xx - (a x H xx + a 2 H 22 + a 3 H 33 ) , (73) 

we get 

3a x H xl - (0-1 Hu + a 2 H 22 + 0-3-^33) 
= (RT) U = #3(01 - 02) + #2(01 - 03) • (74) 

We note however that 

30i-ffii — (01^11 + o 2 H 22 + o 3 H 33 ) = 2a x H xl — o 2 H 22 — a 3 H 33 

= [a x H xx - <J 2 H 22 ] + [a x H xx - 03-^33] 

= [a x H xx + ((7i + 03)^22] + [0i#u + (01 + 02)^33] 

= [a x {H lx + H 22 ) + a 3 H 22 ] + [a x {H xx + H 33 ) + a 2 H 33 ] 

= —aiH 33 + a 3 H 22 — o\ H 22 + o 2 H 33 

= H 33 (a 2 - a x ) + #22(03 - 01) 

= —{RT)\i . (75) 



Thus 



together with (|75|) gives 

= 3a x H xx - (a x H xx 
= {RT) xx = H 3 {a x - 



- a 2 H 22 + a 3 H 33 ) , 
02) + #2(01 - 03) 



(76) 
(77) 



18 



Similar results hold for a = d = 2 and a = d = 3 respectively 

= 3a 2 H 22 - (<riH n + a 2 H 22 + cr 3 H 33 ) , (78) 
= (RT) 22 = H 1 (a 2 -a 3 ) + H 3 (a 2 -a 1 ) . (79) 

= 3^3^33 - (criH u + a 2 H 22 + 0-3^33) , (80) 
= (RT) 33 = H 2 (a 3 -a l ) + H x {a 3 - a 2 ) . (81) 

If we write (|77|),([79]),(PI|) together we get 

= (iZr)ll = #3(<Tl-<7 2 ) + # 2 (<7i-<7 3 ) , (82) 

= {RT) 22 = H l {a 2 -a 3 ) + H 3 {a 2 -a l ) , (83) 
= ( J RT) 3 3 = ^2(a 3 -a 1 )+if 1 (a3-(T 2 ) , (84) 

which are consistent with each other. 

Similarly if we write equations (|7B|) , (|T8"D , (|8"0"|) which are equivalent to 
O,©,© we get 

= 3a 1 H 11 - {a 1 H 11 + a 2 H 22 + a 3 H 33 ) , (85) 
= 3a 2 H 22 - ((TiH n + a 2 H 22 + cr 3 H 33 ) , (86) 
= 3cr 3 # 33 - (0-1 F11 + o- 2 #22 + 0-3-^33) , (87) 

(88) 

from which it follows that 

(7i#ii = G 2 H 22 = a 3 H 33 . (89) 

Equation (p9|), which is the required equation (37), is a condition arising 
from consistency requirement on the u H a d" for a = d. 

Similarly for non-diagonal elements a 7^ b equation (|57|) 

(RT) 23 = \a x H 23 , (RT) 3l = \a 2 H 31 , (RT) 12 = \a 3 H 12 , (90) 

and the time propagation equation (|63| ) yields 

^3 = 0, (a = 2,d = 3); 
a 2 H 13 = , (a = 1, d = 3) ; 

a 3 #i 2 = 0, (a = l,d = 2). (91) 
which is the required equation (34). 
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